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We study the dielectric properties of graphene in the presence of Rashba and intrinsic spin-orbit 
interactions in their most general form, i.e., for arbitrary frequency, wave vector, doping, and spin- 
orbit coupling (SOC) parameters. The main result consists in the derivation of closed analytical 
expressions for the imaginary as well as for the real part of the polarization function. Several limiting 
cases, e.g., the case of purely Rashba or purely intrinsic SOC, and the case of equally large Rashba 
and intrinsic coupling parameters are discussed. In the static limit the asymptotic behavior of the 
screened potential due to charged impurities is derived. In the opposite limit (q — 0, cj — > 0), 
an analytical expression for the plasmon dispersion is obtained and afterwards compared to the 
numerical result. Our result can also be applied to related systems such as bilayer graphene or 
topological insulators. 
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I. INTRODUCTION 

It is now well established that at low energies the 
charge carriers in graphene are de scri bed by a Dirac- 
like equation for massless particles!^ While standard 
graphene, i.e., without any spin-orbit interactions (SOIs), 
does not exhibit a band gap, a gap opens up in 
the spectrum if one includes purely intrinsic spin- 
orbit interactions. 3 The corresponding energy dispersion 
resembles that of a massive relativistic particle with a rest 
energy which is proportional to the spin-orbit coupling 
parameter (SOC). Including SOIs of the Rashba type, 
e.g., by applying an external electric held, lifts the spin 
degeneracy. Depending on the ratio of the intrinsic and 
the Rashba parameters a gap can occur in the spectrum 
or not. 

Many theoretical studies on the dielectric function of 
various systems have been made in the last years. Besides 
semiconductor two-dimensional electron gaseiP^ anc j 
hole gas systems^ large investigations have been made in 
graphene. Starting from the simplest possible graphene 
model within the Dirac-cone approximation^"^ more 
and more extensions have be en in cluded. These 
extensions range from numerica l 11 * 12 ! and analytical 
tight-binding studies and the inclusion of a hnite 
band K a-P to double- and multilayer graphene 

samplesp2H22] graphene antidot lattices,^ and graphene 
under a circularly polarized ac electric fieldP^l 

In this work we study the dielectric properties of 
graphene including both the Rashba and the intrinsic 
spin-orbit coupling. While the case of purely intrinsic 
interactions is well understoodp^KH the dielectric 
function for the general case, where both types of SOIs 
are present, is unknown. Other previous studies have 
investigated the effect o f SOI on magnetotransporPSl and 
the optical conductivity! 2 ** 1221 

Our study is motivated by recent experimental 



and theoretical works demonstrating that the SOC 
parameters can significantly be enlarged by choosing 
proper adatoma^Sni or a suitable environment 



Information that can be extracted from the dielectric 
function range from the screening between charged 
particles to the collective charge excitations formed due 
to the long-ranged Coulomb interaction. Knowledge 
of the latter is not only important for possible future 
applications in the held of plasmonics, where graphene 
seems to be a promising material,^ but also because 
of fundamental reasons. Recent experiments and 
theoretical studies showed that interactions between 
charge carriers and plasmons in graphene, forming so- 
called plasmarons, yield to measurable changes in the 



energy spectrum 
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The paper is organized as follows. In Sec. [TTJ 
we introduce the model Hamiltonian including the 
eigensystem and summarize the formalism of the random 
phase approximation (RPA). In Sec. Ill analytical 



and numerical results for the free polarization function 
of the undoped and the doped system are given. In 
Sec. |IV[ the dielectric function is used to analyze the 
static screening properties due to charged impurities. 
We provide qualitatively the asymptotic behavior of the 
induced potential. The long-wavelength collective charge 
excitations of graphene are derived in Sec. M and 
afterwards compared to the numerical result. We find 
the existence of several new potential plasmon modes 
that are absent without any spin-orbit interactions. Most 
of these zeros, however, are overdamped as can be seen 
from the energy loss function. We close with conclusions 



and outlook in Sec. VI Finally, in Appendixes [A| and |B| 
we give details of the calculation of the free polarization 
function. 



cc 5 





10 

5 






1 1 y 
/ / 


-\\ 


/ / ' 


\\ 


/ / 




/ / ~ 


\ \ 


/ / 


s \ 


/ 


\ 


\. 


/ / 


"~ N \ 


/ / 


\ \ 


/ / 


\ \ 


/ / 


\ \ 


'/ / 


\ \ 


/ 


\ 


Z , 1 , 


1,1,^ 



-10 -5 5 10 -10 -5 5 10 -10 -5 5 10 



klX 



R 



k/X 



R 



k/X 



R 



FIG. 1: Energy dispersion in units of \r for E + ± (solid lines) and E_± (dashed lines): (a) A/ = 2A^, (b) Xj — \r, 
(c) A 7 = X R /2. 



II. THE MODEL 



The eigensystem reads 



We describe graphene with SOI within the Dirac cone 
approximation. At one K point, the Hamiltonian is given 

bj0 



H = v F p • r + X r (t x a) e z + X T r z a z 



(1) 



The Pauli matrices t (a) act on the pseudospin (real 
spin) space. The other K point can be described by 
the above Hamiltonian with a x — ¥ —a x and o z — \ —o z . 
Since the two K points are not coupled, we can limit 
our discussion to the above Hamiltonian, multiplying 
the final results by the valley index g v — 2. Moreover, 
without loss of generality, we assume a positive Rashba 
and intrinsic coupling as the eigensystem and thus the 
dielectric function is not changed for negative values. 



Solution 



For a sufficiently large intrinsic coupling parameter, 
A/ > Xr, the system is in the spin quantum Hall phase 
with a characteristic band gap. For Xu > Xj the gap 
in the spectrum is closed and the system behaves as 
an ordinary scmimctal. At the point where Xr = Xi 
a quantum phase transitions occurs in the system. In 
the following we mainly set vf = 1 and h — 1 . 



\x±±(k)} = 



Ix±t(*0> = 



1 



1 

71 



/ sin(0 T /2) \ 
cos (0 T /2)e tip 
±cos (6> T /2)e lv 
\±sin(6» T /2)e 2 ^y 

/ cos((9 T /2) \_ 
-sin {6 T /2)e lv 
Tsin (6> T /2)e^ 

\±cos (9 T /2)e 2lv ) 



(2) 



with sin ( 
1*1) 



and X a — Xr + aXj, and (k 



E a p{k) — a Xr 



fc 2 + A_ Q 2 (a,0 = ±l). (3) 



For Xr ^ the spin degeneracy is lifted and two distinct 
Fermi wave vectors, k F ± = 

exist. In Fig. [T] the energy dispersion is shown for three 
characteristic values of the SOI. 

The energy scales for the SOC parameters in 
monolayer graphene, Aj = 12/ieV and Xr = 5/xeV for 
an electric field of 1 — , are generally small 36 . However, 
it was shown that these parameters can be enlarged to 
Xi « 30meV for thallium adatomPsl or Xr « 13meV for 
graphene placed on a Ni(lll) surface.^ 

The above Hamiltonian with only Rashba coupling can 
be mapped onto the bilayer Hamiltonian without SOI, 
relating the intcrlaycr hopping parameter tiL ~ 0.2eV 37 
to the Rashba SOC. Our findings can also be applied to 
a topological insulator within the Kane-Mele modelP 
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Dielectric function 



by 



In order to find the dielectric function in RFAPS given 



e{q,w) = 1 - V(q)xo(q,u), 



(4) 



is the Fourier transform of the 

and 



where V{q) = ^ 

Coulomb potential in two dimensions, V{r) — 4 ^ e - 



eo the vacuum permittivity, one needs to calculate the 
free polarization, 



Xofew) = / 

„ LI J 



g v d 2 k 
, 72^)2 

af{E, nm {k)) 



Xr]3V4 



- a [Erisv* ( k + <l)~ E mm ( k )} + i0 ' 



(5) 



In the following we assume zero temperature. The Fermi 
function f(E) then reduces to a simple step function. 
Because of the general relation Xofe~ u ) = [xofe w )]*> 
we restrict our discussions to positive frequencies lu. 



III. ANALYTICAL RESULTS 

A. Zero doping 

For zero doping the valence bands are completely filled 
while the conduction bands are empty. Only transitions 
between bands _E Q _ and Ep + are possible. The resulting 
charge correlation function can be decomposed as 



Xofew) = Y \ 

r)i=± 



VI >»?3 + 



(6) 



Here we introduced the notation -^vim-^vaiu w ) 
describing transitions from the initial band E VlV2 (k) to 
the final band E^ zVi {k + q). For the imaginary part we 
find 



3g 4 



_ 9v 
4A± V - 5g 2 w 2 



4A H 



2lu 4 



(lu 2 - q 2 f /2 
q 2 - uj (u - 2A±)| + \q 2 - lu (oj + 2A±)| 

LU 



and 

lm{ X ± -^ + (q,o J )} = 



[cu 2 ± -q 2 - 4 7 2 ] X 



(7) 



(8) 



u? ± - q 2 - 



-lu(lu±2\-) 
2lu 



lu(lu±2X + )\ 
2cu 



Here we defined u>± = lu ± 2\r and 7 = maxjA^, A/}. 
For equally large spin-orbit coupling parameters, A^ = 
A/, the imaginary part is divergent at the threshold 
lu = q but finite otherwise. The divergent part of the 
polarization is x H >++ as the bands E + ±(k) are linear 
in momentum. 

The real part can be obtained via the Kramers-Kronig 
relations 



Re{xo(g,w)} = -V I dcu . (9) 



After carrying out the remaining integration, where it is 
necessary to keep the principal value, we arrive at 
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FIG. 3: Static real part of the charge susceptibility for (a) undoped graphene with fixed A+ for Xr — 2Xi (dotted), 
A_r = Xi (dot-dashed), A/ = 2Xr (dashed), Xr = Xj = (straight), and doped graphene with (b) Xr = Xi = 0.3/x, 
(c) Xr = 2Xi = 0.3/i, and (d) 2Xr = Aj = 0.3/1 in units of the density of states D(0) = g v /z/7r. 
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FIG. 4: Asymptotic screened potential (in units of QA+/eo) of undoped graphene for fixed A + = A/? + A/ and 
different spin-orbit coupling parameters, (a) \r — A/ (straight line), A/ = 2 A/? (dashed), Ajj = (dotted). Also 
shown the non-interacting case A^ = Xj = (dot-dashed), (b) \j = (straight), A^ = 2A/ (dashed). 
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Here we defined the function 

22 22 
jCa (^J = x + A In s In 



2lo 



X + U! 



(12) 



is analogous). The free polarization in the doped case 
reads 



B. Finite doping 



Xo(g,w) = xo(q,oj) + 8xk F+ (q,u) +(5xfc F _(g,w). (13) 



We now continue with the case of a finite chemical 
potential lying in the conduction band (the p-doped case 



Xo is the undoped part given above. The two remaining 
contributions, $Xk F+ and Sxk F ^ , with 



$Xk F± (q,u) 



9v \ " 

4tt 2 ^ 

a,fi,iy—± 1 



v 



k F ± 



d 2 k £ 

a = ± 
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(x±+(k) 


X^uik + q)) 


2 


oj + iO — a [E^ v 


{k + q) - E ±+ (k)} 



(14) 



refer to transitions with initial states in band E ++ and 
respectively. As the expressions for the extrinsic 
real and imaginary part of the free polarization function 
are quite lengthy, we refer to Appendix [B] where the 
results including major steps of the derivation can be 
found. 



Similar to the undoped case, the density correlation 
function of graphene is finite at uj = q for A# 7^ A/ and 
divergent for A/? = A/. However, this divergence vanishes 
in the RPA improved result P 

From the shape of the Fermi surface and the dispersion 
relation as given in Eq. (§, we can determine the 
boundaries of the dissipative electron-hole continuum.^ 



G 



a) 



2x10 
1x10 



-1x10 
-2x10 



c) 



2x10 
1x10 



-1x10 
-2x10 




b) 


2xl0" 2 




lxlO" 2 




CM 




lxlO" 2 




2xl0" 2 


d) 


-2 

2X10 


x> 


-2 

IXlU 








3 


lxlO" 2 



-2x10" 



"i — 1 — r 



-i 1 r 




J i L 



20 



40 
lir 



60 80 




FIG. 5: Asymptotic screened potential (in units of Q/i/eo) of doped graphene for various spin- 
parameters: (a) (Afl/jti, A///x) = (0,0.3), (b) (0.3,0.15), (c) (0.15,0.3), and (d) (0.3,0.3). 



In Fig. [2] this is shown for the particular choice of 
Xb = 2Xi = 0.3/x. In general, the lower and upper 
boundaries of the damped region I are given by 



^low = max |0j v ( k F - -q) +x\- x /k 

and 



2 1 \2 
F- + A + 



US 



z up = max I \J{k F _ +qf + X\ - yj k 2 F _ + X 2 + , 
yj (k F+ + qf + XI- y/k% + + XL ] j , 



respectively. Region I is due to intraband transitions 
from band E± + (k) to E±+(k + q). Region // accounts 
for interband transitions between conduction bands and 
is confined by 



Up/ LC 



= d(k F - ± q y + x 2 _- Jk F _ + x 2 + + 2x R 



For region III the lower limit reads 



low 



(fc F _ - qf + X\ + Jk 2 F _ +X 2 _- 2X R , 



{k F+ - q y + x 2 _ + Jk 2 F+ + x 



while there is no restriction to the upper boundary. This 
part is due to transitions between valence and conduction 
bands. 



IV. SCREENING OF IMPURITIES 

The potential of a screened charged impurity is 
obtained from the definition of the dielectric function, 



$(r) = 



Q 



£ o Jo 



dq 



Jo(q r ) 

e(q,oy 



(15) 



Jo (a;) is the Bessel function of the first kind and Q the 
charge of the impurity. Making use of Eq. fi"5"l ) the 
screened potential for the undoped system is calculated 
numerically where is mainly determined by the long- 
wavelength behavior of the static correlator!^ As can be 
seen from Fig. ^a), the long- wavelength limit of the 
polarization xo(0, 0) is finite in the semimetallic state 
(Xr > Xi) and zero otherwise while for large momenta 
all functions scale like 1/q. From Fig. |4|a) we can see 
that for Xj > Xr the potential scales like $(r) oc 1/r at 
large distances. For Xr > A/ the asymptotic potential 
behaves as $(r) oc 1/r 3 ; see Fig. |4^b). The actual values 
of fir at which the above asymptotics are appropriate 
approximations depend on the difference of Xr and A/. 
As mentioned in the introduction, the two different 
parameter regimes belong to different phases separated 
by the quantum critical point at Xr = Xj. 

The static density correlator for the doped system 
is much more complicated. Integrals of the form ( 15 ) 



are usually treated analytically by approximating the 
Bessel function by its asymptotic values. The subsequent 
Fourier integral can then be solved with the Lighthill 
theoremPS The above theorem states that singularities 
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FIG. 6: Density plot of Im{— l/e(q,uj)} for various spin-obit coupling parameters: (a) (A«//i, Xi/fi) = (0,0), (b) 
(0,0.5), (c) (0.25,0.25), (d) (0.5,0). Straight lines correspond to the numerically calculated zeros of the real part of 
the dielectric function while dashed lines represent the long- wavelength result of Eq. (19 1. 



in the derivatives of the dielectric function give rise 
to a characteristic, algebraic, oscillating decay of the 
screened potential. Physically, these Friedel oscillations 
are due to backscattering on the Fermi surface. We 
can thus make qualitative predictions for the potential 
<I>(r) at large distances away from the impurity, only 
from the analytical structure of the polarization function 
without carrying out the integration. Afterwards these 
predictions are compared to the exact numerical solution. 

For non zero SOC and Xr ^ A/ the first derivative of 
the polarization function is singular at special points q = 
2kp±] see Figs. [3]jc) and (d). According to the Lighthill 
theorem the potential will exhibit a superposition of 
two different kinds of oscillations whereat $(r) oc 1/r 2 . 
This beating should be observable in sufficiently clean 
samples if the Rashba parameter, and the consequential 
breaking of the spin-degeneracy, is large enough. For 
predominant intrinsic SOI, the two oscillatory parts 
interfere constructively finally yielding an additional 
spin-degeneracy factor of g s = 2™ For A# = A/ already 



the first derivative of XofeO) is singular at q = 2kp- 
while at q = 2kF+ only the second derivative diverges; see 
Fig. |3jb). The main contribution in the potential again 
will be of order 1/r 2 . The numerical inspection of $(r) 
as displayed in Fig. [5] confirms the above predictions. 

The resulting potential deviates significantly from the 



$(r) cx 



cos (2fci?r) 
(2k F rf 



(16) 



behavior of st anda rd graphene within the Dirac cone 
approximat ion jSEHl 

Nevertheless, including the full dispersion of graphene 
can also lead to a different decay behavior, i.e., to 
anisotropic regular Friedel oscillations decaying like 
l/r 2 P 
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FIG. 7: Energy loss function Im{— l/e(q,uj + iO)} for fixed q = O.lfi with (a) A_r = 0, Xi = 0.5/x and (b) A^ = 0.5^, 
A/=0. 



V. PLASMONS 

Plasmons are defined as the zeros of the dielectric 
function, 



e(q,u p - ii) = 0. 



(17) 



For small damping constant 7, Eq. (17 1 can be 
substituted by the approximate equation^! 



Re{e(q,uj p )} = 0. 
Only if 7 is small compared to 



(18) 

one can speak 
of collective density fluctuations. For large Landau- 
damping, it is thus important to also discuss the more 
general energy loss function Im{— l/e(q,uj)} which gives 
the spectral density of the internal excitations of the 
system. „„ 

Similar to Refs.LLiiLS there are several solutions of Eq. 
( 18 1 for non zero SOC parameters. In Fig. |6j these 



solutions are shown as straight lines together with a 
density plot of the energy loss function Im{— l/e(q, ui)}. 
One of these solutions has an almost linear dispersion 
with a sound velocity close to the Fermi velocity which 
exhibits an ending point for Xr ~ A/ associated with 
a double zero of the real part of the dielectric function. 
However, as can be seen from Figs. [7]ja) and (b), this 
solution does not yield to a resonance in the loss function 
and does thus not resemble a plasmonic mode. In the 
case where the gap in the spectrum is closed (Xr > A/), 
two additional zeros appear leading to potential high 
energy modes similar to bilayer graphene! 18 * 19 ! However, 
these potential collective modes are damped by interband 
transitions; i.e., the corresponding peaks in the loss 
function are broadened out as can be seen from Fig. gb) 
and no clear signature is seen in the density plot. 

We are thus left with the branch which is also present 
for "clean" graphene and which resembles the only 



genuine plasmonic mode; see Fig. [6|a). Its dispersion 
ujp can be approximated in the long-wavelength limit 
(q < u) by32l 



u°(q) = 0Jq, 



(19) 



where the prefactor is given by /3 = 

We thus recover the typical 



3f "fy 

Y^-dispersion of 2D-plasmons. 

The long-wavelength approximation is shown as a 
dashed line in Fig. [6] and coincides with the numerical 
solution, uj p , for small momenta, whereas for larger 
momenta, the ui p is red shifted compared to u)®. If 
Xi is large enough, u p remains in the region where 
Landau damping is absent, see Fig. j^b),™ otherwise 
it eventually enters the Landau-damped region due to 
interband transitions from the valence to the conduction 
band, see Figs. [6]ja), (d). 

For two occupied conduction bands, which is the case 
in Fig. |6jc) and in Fig. [8j the plasmon mode is disrupted 
at q ss 0.05/z by a region with a finite imaginary part 
where it becomes damped. This additional Landau- 
damped region is due to interband transitions from the 
two conduction bands. The analytical description of the 
boundaries of this region can be found in Sec. |III|B. 



This "pseudo gap" of the plasmon dispersion can also 
be obtained from only considering Eq. (18) since the 



"plasmon" velocity formally diverges at the entering and 
exit point as can be seen from Fig. [6](c). The crossing 
points can alternatively be approximated by looking 
at the intersection of this region with the analytical 
long-wavelength approximation of the analytic plasmon 
dispersion. For the quantum critical point (Xr = A/), 
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FIG. 8: Energy loss function Im {-l/e(q, w + iO)} for (a) (A H //z, Aj//z) = (0.25,0), (b) (0.25,0.25), (c) (0.25,0.5), 
and (d) (0.25,0.75). The straight blue lines show the undamped plasmon modes. The black lines indicate the 
boundaries of the single-particle continuum (see Sec. Ill B). 



this leads to the critical wave vector 



(p - \jp 2 T 4 (kZ - yfk 2 F _ + 4Af + 2A fl ) 



and in particular to q~ r w 0.019/i and g+ « 0.025/i 
for A fl /j = 0.25/1. For a proper analysis, the full 
energy loss function thus needs to be discussed which 
is done in Fig. [9| It shows how the spectral weight is 
eventually transferred from the lower to the upper band 
as momentum is increased, explaining the step in the 
plasmon spectrum as shown in Fig. [Hj 

The pseudo gap of the plasmonic mode always appears 
for Xfj < 0.5/1, since then two conduction bands 
are occupied independently of the value of A/, but it 
decreases for increasing A/ as the dissipative region due 
to interband transitions diminishes. In the opposite case 
of Xji > 0.5/x, either one or two bands can be occupied. 
For zero intrinsic coupling, the pseudo-gap is absent but 
increases up to a maximum value at around Xj rs A^ for 
increasing A/. 

Let us close with a comment on plasmons in undoped 
graphene. For neutral monolayer graphene and at zero 



temperature, plasmons can exist if one takes into account 
a circularly polarized light fielcP^ or effects beyond 
RPAp^ and in bilayer by including trigonal warping. — 
In our system, the real part of the dielectric function 
is always nonzero for the undoped case and thus no 
plasmons exist. 



VI. CONCLUSIONS AND OUTLOOK 

We have presented analytical and numerical results 
for the dielectric function of monolayer graphene in the 
presence of Rashba and intrinsic spin-orbit interactions 
within the random phase approximation for finite 
frequency, wave vector, and doping. The cases of 
predominant Rashba and intrinsic coupling and the case 
of equally large SOC were opposed. 

In the static limit the screening properties due to 
external impurities were studied. Our findings show 
that the power-law dependence of the screened potential 
in the undoped system depends on the ratio of the 
Rashba and intrinsic parameters. While for A# > Xj the 
screened potential scales like <&(r) oc 1/r 3 , for X] > Xr 
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FIG. 9: Top: Energy loss function Im{— l/s{q, w + iO)} for A# = 0.25/x = Xj = 0.25/1 and various wave vectors q. 
Bottom: The same for \r = 0.25// and A/ = 0. 



a weaker screening, $(r) oc 1/r, was found. For finite 
Rashba coupling, a beating of Friedel oscillations in the 
doped system occurs due to the existence of two distinct 
kinds of Fermi wave vectors. For large A/ 3> A^, 
this beating vanishes and the two contributions interfere 
constructively. 

In the last section the influence of SOI on the 
collective charge excitations was discussed. We found 
that while only one plasmon mode exists for standard 
graphene, several new potential modes occur for finite 
SOC. However, most of these modes are overdamped 
and can hardly be detected as they lie in the region 
with finite Landau damping. In the case when the two 
conduction bands are filled, the undamped plasmon mode 
is disrupted by a narrow dissipative region strip due 
to particle-hole excitations. This "pseudo gap" might 
be useful to gain further control in possible plasmonic 
circuitries based on graphene. 

As already mentioned in the beginning, our findings 
go even beyond monolayer graphene. For purely Rashba 
coupling the dielectric function presented in this work 
equals that of bilayer graphene. The role of the SOC 



parameter is then played by the interlayer hopping 
amplitude tn, being several orders of magnitude larger 
than Afl. Additionally, the Hamiltonian in Eq. ([I]) 
generally describes a system known as the Kane-Mele 
topological insulator!^ Our discussion can thus be fully 
adopted to materials modeled by this Hamiltonian. 
Besides that, our findings might also be relevant for other 
monolayers with similar symmetry properties compared 
to those of graphene, e.g. MoS 2 , where SOC is naturally 
strong™ A detailed study of the dielectric properties of 
M0S2, however, is left open for future works. 



Acknowledgments 

We thank G. Gomez-Santos and O. V. Gamayun 
for useful discussions and comments. This work 
was supported by Deutsche Forschungsgemeinschaft via 
Grant No. GRK 1570, by FCT under Grant No. 
PTDC/FIS/101434/2008, and MIC under Grant No. 
FIS2010-21883-C02-02. 



11 



Appendix A: Details of the calculation of the undoped polarization 



The undoped polarization is composed of four parts, 



(Al) 



As two of them can be obtained by a simple substitution, i.e., X\ R ~^ \q,uj) — X—x~* + (1t uj ) an d X\ a ~^ (q,uj) = 
XZ\^ ^(q^) only two contributions remain. With the help of the Dirac identity the imaginary parts read 

Imjx >- + (<Z,w)} = -f L [d 2 k aU X -^(k)\x-+(k + q))\ 2 S[u J -a(E^ + (k + q)-E^(k))} (A2) 

** Li 

(A3) 



a=±l 



g^ 2 -, 2 -4A + 2 ] 



3q 4 -4A+V - 5q 2 uj 2 + 2uj 4 \q 2 - u (u) - 2A+)| + \q 2 - uj (u) + 2A+) | 



(uj 2 - g 2 ) 3/2 



and 



Im{ X + -^- + (q,u)} = -^- f d 2 k ]T a\(x+-(k)\x-+(k + q))\ 2 5 [uj - a (E_ + (k + q) - E + _(k))} (A4) 

•* „ LI 



a=±l 



8tt 



dy 



2 - / g 2 (g»-^+4Af)( g '-^+4A' H ) 



|A_| 



(y- A_)(w-y + A_) 



1 - 



cj 2 - q 2 - 2uj + y - 4A J? Aj 



9v 



e [co 2 ± -q 2 -±i 2 ] 



e [uj% -q 2 - 4 7 2 ] 

2qJy 2 -X 2 _ 

q 2 -cj(w±2A_)| U 2 -w(w±2A + ) 



c4 - <? 2 



2w 



2uj 



with 7 = max{A^, A/} and y = w fc 2 + A 2 . 

We can now make use of Eq. ([9| in order to find the real part. The first contribution reads 

2 f°° uj' 

Re{ X >- + (<7,^)} = -V \ diu'— -Im{ X ^(q^')} 

k J i J uj' z — UJ A k J 

l - {-A A (4A 2 ) + £ A (vV+4A 2 ) + 2£_ A ( V /^TA 2 ~+ A) - £_ A (\Az 2 + 4A 2 ) 

Zq 2 u? + 4g 2 A 2 - 3<? 4 - 2w 4 tt ' 



8tt 



+0 [q - uj 
Here we introduced the functions 



(q 2 - UJ 2 ) 3 / 2 



K\{x) = 2^fx + 



2 A 2 



arctan 



(q 2 - UJ 2 )yJx 



- (3g 4 - bq 2 uj 2 + 2uj 4 - Aq 2 X 2 ) Re 



(q 2 - w 2 ) 3 / 2 



(A5) 
(A6) 



(A7) 



(A8) 



(A9) 



and 



C\ (x) = x + A In a; — uj — 



a ,. 2 | - 2 -? 2 ln 



2w 



X + UJ 



X — UJ 



(A10) 
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The second contribution can be solved in a similar way, 

Re{ X + -^- + (g,w)} = -V \ dJ— -Im{ X +— + (q,oj')} 



9y_ 
47T 



2X R (1 + In 4) - - Re I ^q 2 - (lo + 2X R ) 2 arcsin 



uj + 2X r 



\J q 2 — {—>jJ + 2X R ) 2 arcsin — ^ | 



-- [G f+Xa (vV+4 7 2 - w - 2A«) + G_« +a r (vV + 4 7 2 + w - 2A«)' 
+0 [A fl + A/] £ A _ ^ + A 2 _ - A -) - ^ si g n + A /) ^- (vV+4 7 2 - 2A fi ) 
+0 [A fl - A,] £ A+ (Vg 2 + X\ «gn (A« - Aj) £ A+ (vV + 4 7 2 - 2A«) 



(All) 



with 



G w (x) =y / (a; + w) 2 - q 2 + uj In (y^ + uj) 2 - q 2 + x + uj) 

/— ujx + uj 2 - q 2 + yJuS 2 - q 2 J (x + uj) 2 - q 2 

- v uj 2 - q z In 



(A12) 



Appendix B: Details of the calculation of the doped polarization 



The extrinsic part for the band E h , 



5xk F Aq^) = £ 2 



'- r ? r 



. } „ - '\X-+(k) X^{k + q) 
d k J^ 1 oj + iO-a[E llv (k + q)-E_ + (k)Y 



(Bl) 



can be summarized as 



<5x/c F _ (q, w) 



= 9v_ v f kF 
^ 2 Jo 



d 2 k 



+ zO+ v /fc2 + A 2 h + ^|fc + q| 2 + A 2 h ^ |( X _ + (fe)| X — (k + q) 
(u + iO + ^k 2 + \ 2 + y - (jfc + q| 2 + A^) 



+ 



+ 



+ 



(u + iO + ^k 2 + X\ - y/\k + q\ 2 + X 2 ^j |( X -+(fe)|x— (k + q) 
(oj + tO+^k 2 + X 2 + y - (\k + q\ 2 + X 2 + ^j 

(w- + *0 + ^k 2 + X\ - yftk + q\ 2 + X 2 ^ |(x-+(fc)|x+-(fe + Q) 
(w_ +z0 + ^fc 2 + A2.) - (jfc + gl' + A 2 ^ 

(w- + «0 + ^fc 2 + A 2 , + + q| 2 + A 2 _) |(x-+(fe)|x++(fc + q) 



+ (w -> -w) 



2 

(w_ + iO + ^k 2 + \%) - (Jfc + q| 2 + A 2 ^ 



(B2) 



where (w — > — uj), and thus (w_ — > — denotes terms with the sign of the frequency changed compared to the 
preceding expression. The corresponding expression for E ++ can be obtained by substituting X R — > — A# and fcp_ — > 

kp + . After carrying out the angle integration for the real part and choosing a proper substitution, x = , Ik 2 + X 2 , — A + , 
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we arrive at 



Re{* X fc,-(?,w)} = -g : R(^-p /" 



x + X R 
2x 



+ 



[q 2 -jx + % ){x + | + A + )] 2 sign (q 2 - u? - 2lo{x + A + )) 



x{x + %) 



4A 5 ; 



4x(x + us) 

ign (q 2 - <jj 2 _ - 2cu-(x + A+) - 4A fl A/)] + (w -w) | 



x si 



(B3) 



These integrals can now be solved in terms of trigonometric and hyperbolic functions.^ In order to simplify the 
expressions we use the shorthand notation^ 



f(x) = sign (b - x) (f(b) - f{x)) - sign (a ~ x) (f(a) - f(x)) 



(B4) 



The result can then be written as 

Re {S XkF _ ( q , u )} = _ ln M - A, 



27T 4-7T 



" A ' /V +u; 2 - 2.x A, 



2u 



-sign(u;_) [H% 



q 2 -u 2 + 2A_w 

2777 



g 2 +cj 2 - 2A+w 
2777 



fi — Xi -\-cu 



+ (u> -» -u>) (B5) 



with 



a" 2a? + Bfx + ZJoffJr? g? / , , \ 

, ln v . v , + ln ^ 2 y 7r + ^ + (B6) 



and = a" + Bfx + jfx . The coefficients read 

, 2 



{q 2 -lo{lo + 2X + )Y 



> "2 



{q 2 - w (w_ - 2 sign (w) A/))' 



# = (w 2 - q 2 ) (2A + + u) , /3£ = 8\ 2 R \i -2\ R \ I (u> + \u\) - sign (u) {to - 2\ R ) (q 2 ~ ui 2 _) 
a a 2 u 2 a _ x 2 2A + + 7c 4A^ - 4g 2 + 8A^ + 5c 2 2< 

7x=«-ff , 7 2 =M_-<? , c i=y + ^2— + § w 



a 



? (2A+ - u) 3lo 3 + 6X 2 + uj 2 - 4(q 2 + \ 2 + )lu - 8X 3 + 



72 



4 7 ^ 



16 



) L 2 



( <-2 



The calculation of the imaginary part is quite similar. Starting from Eq. (B2 1 and carrying out the angle integration 
in a way similar to the real part, we arrive at 



Im{6 XkF M")} = ^ R& { V 



H—Xi 



dx 



[ql-(x+%){x+± + \ + )X 



sign (x — A_ + us) 



x(x+ |) 



3ZI1+ 4A + 
" 1 1 + q 2 -u 2 



4a;(x + w) 



sign (ac + A + + w) 



-( 2 ;+|+A + )2 

(w -» — w) I 

(B7) 
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The result can again be written as 



Aj+LJ 



-u] 9 \/i -Xi + oj] sign (A+) [B% (-e) - i^(e)] + iK% (A_ - w) 



^2 (A-) 



l-L — Xj +CJ 



-i0[-u;]0[^-Aj + u;]sign(-A_)[E£(-e)-i?£(e)]} - (w -> -u>) 



(B8) 



The limit e — > in Eqs. (|B5j) and (B8) can now be taken safely giving finite results. 
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address: andreas.scholz@physik.uni-regensburg.de 
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